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We demonstrate a proof-of-principle magnetometer that relies on the active oscillation of a cold atom Raman 
laser to continuously map a field-sensitive atomic phase onto the phase of the radiated light. We demonstrate 
wideband sensitivity during continuous active oscillation, as well as narrowband sensitivity in passive Ramsey- 
like mode with translation of the narrowband detection in frequency using spin-echo techniques. The sensor 
operates with a sensitivity of 190 pT/-\/Hz at 1 kHz and effective sensing volume of 2x 10~'^ mm'^. Fundamental 
quantum limits on the magnetic field sensitivity of an ideal detector are also considered. 



Magnetometers are useful in a broad range of appli- 
cations, including NMrf^, fundamental physicd^, bio- 
sensinj^, and atmospheric physic^. A wide array of 
physical systems with various trade-offs in sensitivity, 
bandwidth, size, and operating conditions have been em- 
ployed. Examples include magneto-resistive materials^, 
squids'^, and nitrogen- vacancy centers in diamond^. 
Atom-based sensors offer several advantages: intrinsi- 
cally high sensitivity to magnetic fields, well-determined 
absolute calibration and scale factors, no need for a cryo- 
genic environment, and remote sensing capabilities via 
optical probinglSHlll 

In this Letter, we demonstrate a unique atomic magne- 
tometer based on Raman lasing transitions between hy- 
perfine ground states of an ensemble of 10^ laser-cooled 
^''Rb atoms confined in a low-finesse optical cavitjff^lti^. 
The lasing frequency is sensitive to magnetic fields rather 
than the cavity resonance frequency because the laser 
operates deep in the bad-cavity, or superradiant, regime 
where the atomic coherence damping rate "fj_ is much 
smaller than the cavity damping rate k. 

The vast majority of atomic magnetometers sense an 
atomic response to magnetic fields through the polariza- 
tion rotation or phase shift of probe light passed through 
an atomic vapoi'^^'^^Hlo] contrast, here we directly 
determine the phase response (f){t) of an atomic dipole 
to external magnetic fields from the phase of the optical 
radiation ip{t) emitted by the dipole. 

We demonstrate that our sensor can transition be- 
tween periods of active oscillation and passive Ramsey- 
like phase evolution, enabling the detection of magnetic 
fields free from possible perturbations, and with high rep- 
etition rates as was first discussed in Ref. ^16, This flexi- 
bility is realized by controlling the optical radiation rate 
via the intensity of a Raman dressing laser. Unlike typ- 
ical good-cavity lasers, the passive mode of oscillation is 
possible in this bad-cavity laser because the atomic gain 
medium is the primary reservoir of phase information. 

In the following, we will also consider two fundamen- 
tal limitations on the magnetic field sensitivity of an 
ideal superradiant magnetometer: quantum Schawlow- 
Townes phase diffusion of the atomic phase (f>{t) at low 
frequencies, and photon shot noise (PSN) associated 
with dete cting the emitted light phase ip{t) at higher 
frequencie d^^l^^l We demonstrate that operation in the 
passive mode, combined with spin echo techniques, allows 




FIG. 1. (a) Experimental schematic and (b) energy level dia- 
gram. The laser-cooled atoms are trapped between the cavity 
mirrors, and the coils are used to apply bias and test mag- 
netic fields. Atoms are continuously optically pumped back 
to It) using two lasers applied from the side (green, violet). 
The dressing laser (red) induces a Raman decay to IJ,), with 
collectively enhanced emission into the cavity mode (blue). 
Both the trapping optical lattice and dressing beams are fil- 
tered from the cavity output before the light is detected at the 
photodiode using a heterodyne reference beam derived from 
the same laser as the dressing laser. 



the enhanced sensitivity of the low frequency detection 
band to be translated to higher frequencies. In principle, 
operation in the passive mode also avoids noise limits set 
by fundamental quantum phase diffusion, and may the- 
oretically achieve sensitivity at the standard quantum 
limit on phase estimation using coherent spin states^^. 

The superradiant magnetometer described here oper- 
ates in a similar manner to the bad-cavity laser previ- 
ously used to explore physicJ^^^^ underlying proposed 
^ 1 mHz linewidth laserg2i"22]^ ^ simplified experimen- 
tal diagram and energy level diagram is shown in Fig. [T] 
The atomic gain medium of the laser consists of = 1 to 
2 X 10^ ^^Rb atoms laser cooled to 40 /iK and trapped in 
a one-dimensional optical lattice at wavelength 823 nm. 
The laser cavity has a finesse F — 710 with mirror sepa- 
ration L = 1.9 cm and TEMqo mode waist Wo ~ 71 /im. 
The root-mean-squared (rms) extent of the cloud is ap- 
proximately 1.5 mm along the cavity axis and 15 ^m per- 
pendicular to the cavity axis, yielding an effective field 
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sensing volume of 2.1 x 10 ^ mm'^. 

Previous studies of steady state superradianc^i^lt^^ 
have utilized the pseudospin-1/2 system composed of 
the magnetic-field-insensitive clock states. Here, we 
use the first order magnetic field-sensitive states |t) = 
|525i/2,i^ = 2,mF = 2) and |;) = \5^Si/2,F = l,mf = 
1) to form the pseudospin-1/2 system. The atomic dipole 
of this two-level system only very weakly radiates at 
the microwave hyperfine frequency difference between 
the two states uthi/^iT = 6.834 GHz. To induce an ef- 
fective optical decay from | f) to 1 1) at a controllable 
single-particle rate 7, we apply a 7r-polarized Raman 
dressing laser beam that is non-resonant with the op- 
tical cavity and detuned A/2tt = 1.1 GHz from the |t) 
to |i) = |5^Pi/2,-F' = 2,m'p = 2) transition with wave- 
length 795 nm. The cooperativity parameter of cavity 
QED, or Purcell factor, is C 7.7 x 10"^ foj. ^i^^ ^^•^ ^o 
14,) Raman transition. 

A TEMqo mode of the optical cavity is tuned within 5 
MHz of resonance with the emitted Raman photon fre- 
quency. Single particle scattering into the cavity is too 
weak for useful sensitivity, with a total scattering rate 
from the ensemble scaling as NCj. However, when the 
total single particle scattering rate into the cavity exceeds 
the atomic coherence dephasing rate j± , the atoms spon- 
taneously synchronize their dipole moments, leading to 
a stimulated enhancement of the emission rate scaling as 
N'^C'y as in the phenomenon of superradiancd^. The 
phase of the collectively emitted light ip{t) adiabatically 
follows the collective dipole moment of the ensemble <j){t) 
within the bandwidth k/2, where k = 27r x 11.1 MHz is 
the cavity power decay rate. Given the resonance con- 
dition and dipole selection rules, the only available state 
for decay through emission into the cavity mode is the 
14,) state such that there is no gain competition between 
different decay paths. 

To operate in a quasi-continuous mode, the atoms are 
continuously recycled at single-particle rate W from \\.) 
to It) by applying two cr^-polarized optical repumping 
beams at 780 nm such that the only dark state is as 
shown in Fig. [TJd. 

The phase noise of the 60 kHz FWHM linewidth Ra- 
man dressing laser is removed from the signal by detect- 
ing the emitted light in heterodyne using a 40 fiW beam 
also derived from the dressing laser. The signal is quadra- 
ture demodulated and both demodulation channels are 
recorded in the data acquisition system to reconstruct 
the full light phasor £{t) = A(t)e*'^(*). 

A set of Helmholtz coils applies a dc bias magnetic 
field Bq — 2.4 X lO^"' T in the ^-direction perpendicular 
to the cavity axis. This establishes the quantization axis 
and shifts the transition frequency by uj^d^'K = 5.1 MHz 
relative to the zero-field ground state hyperfine splitting. 
The dc field also breaks the degeneracy of the mp lev- 
els so that microwave rotations can be performed exclu- 
sively on the It) and |t) states with characteristic Rabi 
frequency J7/27r 35 kHz for the spin echo studies to 
follow. 
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FIG. 2. (a) The power spectral density (PSD) of E" (t) in the 
presence of modulation at 6.9 and 10.2 kHz, averaged over 
many trials. Dashed line shows the expected PSD from the 
applied modulation. The inset shows the detected (points) 
and predicted (line) light phase i/)(f) versus time, (b) The 
calculated magnetic field sensitivity (red) versus frequency 
from detection of 7/)(t) during active laser oscillation. Also 
shown are the expected Schawlow-Townes phase diffusion 
limit (green), the (PSN) limit (black) for g = 7.2 x 10"^, 
and the PSN Umit for <? = 1 (blue). 



We investigate the system's response to small magnetic 
fields by applying small modulations B(€) to the current 
through the same coils, with phase, frequency and am- 
plitude set by the data acquisition computer that also 
records the generated light signals. The atomic transi- 
tion energy E{€) and hence the time rate of change of 
the collective atomic phase depends on the small modu- 
lation of the field through 



E{t)/h-- 



d(j){t) 

dt 



LUM + 2Tra(Bo + B{t)) (1) 



where a ~ 2.1 x 10^*^ Hz/T, and h is the reduced Planck 
constant. 

We apply a small sinusoidal modulation of the field 
B{t) — B,n cos{2Tr f„it + 9m) to create a modulation of 
the atomic phase as (l){t) = /3 sm{2TTfmt + 0„i), where the 
modulation index is /3 = aBm/ fm- Here and in the rest 
of the paper, we ignore the contributions to 4>{t) from the 
hyperfine splitting and the dc bias field. 

Fig. [2] demonstrates the simultaneous detection during 
continuous superradiance of two applied discrete mod- 
ulations at frequencies 6.9 and 10.2 kHz. Also, we ob- 
serve no modulation at the difference or sum of the two 
modulation frequencies, confirming the linearity of this 
detection scheme at the current level of sensitivity. 
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The data was obtained by detecting the Hght phasor 
£{t) during a time window of 3 ms. Atom loss causes 
the superradiance to cease after approximately 20 to 
100 m^^, after which the atoms are reloaded and the 
experiment is repeated every 1 second. The observed 
carrier frequency chirps as atoms are lost due to cavity 
pulling of the oscillation frequency combined with atom 
dependent shifts of the cavity mode frequencjl^^. To 
avoid spectral components from the chirping, we apply a 
third-order digital Bessel high pass filter to the recorded 
phase ip{t) with pass-band cutoff frequency 300 Hz to 
produce a filtered time phase ip'(t). To remove any pos- 
sible amplitude modulation (AM) sideband noise, we set 
A{t) = 1 in the data analysis to construct £'{t) = e^^ 
This could be accomplished in real time by phase lock- 
ing an oscillator with low AM noise to the emitted light 
phase. Spectral leakage of the carrier due to a square 
sampling time window is mitigated by multiplying the 
phasor with a fourth-order Blackman-Harris time win- 
dow wit) such that E"{t) = w{t)£'{t). 

The double sideband power spectral density of phase 
fluctuations S^{f)— is calculated from the Fourier trans- 
form of £"{t) and averaged over many trials. The spec- 
trum of the equivalent rms magnetic field noise density, 
or sensitivity, ri{f) (units of Tesla/VUz) is then calcu- 
lated as 

V{f)^-JSM)- (2) 

Fig. [2]3 shows the calculated ideal detection sensitivity 
and the measured field sensitivity for 7 = 37 sec^^, 
quantum efficiency q — 7.2 x 10^'^, = l.f x 10^, and 
W = 4.5 X 10* sec~^. In our unshielded apparatus, am- 
bient magnetic field noise, noise in the coil driver elec- 
tronics, and emission frequency chirping contribute noise 
far above the Schawlow-Townes limit at frequencies < 11 
kHz. At high frequencies, PSN limits our ability to re- 
solve ac magnetic fields. 

In the absence of technical noise, the ideal phase noise 
density S^^{f) = C"f/{nf)'^ + l/{2md) is expected to 
be fundamentally limited at low frequencies by quan- 
tum phase diffusion of the collective Bloch vector (the 
Schawlow-Townes limit, first term) and at high frequen- 
cies by PSN (the second term). For photon quan- 
tum detection efficiency q and optimum repumping rate 
I^pk — NCjfP^, the rate of detected photons is rhd — 
qRN'^C'j/8. Here, i? is a reduction factor for the multi- 
level ®''Rb scheme presented with value R < 3/5. The 
ideal field sensitivity can be written as: 

where the corner frequency is given by /o = 
y/qRNC'y / {2tt). The minimum detectable field scales 
with the single particle decay rate into the cavity C7, 
while the corner frequency scales with the collectively 
enhanced single particle scattering rate into the cavity 
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FIG. 3. (a) (Top) Time traces of the measured (points) and 
fitted (solid line) B{t) and (bottom) the measured (solid line) 
light phase ^p{t) and expected atomic phase (t>(t) (dashed black 
line) in the presence of magnetic field modulation. During the 
dark period in which all optical pumping and dressing lasers 
are shut off (light blue region), spin echo vr-pulses cause the 
unobserved atomic phase (j){t) to coherently increase in re- 
sponse to the applied magnetic field modulation. This mani- 
fests as a discontinuous advance in the phase of the light when 
the optical pumping and dressing lasers are turned back on. 
(b) The magnitude of the phase advance Ai/j — 4'(Td) — V'(O) 
(points) increases with the number of vr-pulses n-^. The total 
dark time is Td — Ut^t. Also shown is the expected slope of 
2/3 (line). 



NC^. In principle, the tunable decay rate 7 can be re- 
duced until the single-particle transition broadening de- 
scribed by a transverse coherence decay rate 7_l is no 
longer negligible compared to setting a minimum 

C7 ~ for which {if{f)f ^ {l^/N){l + PIP,) 

with corner frequency /o ^ 7_l. These ideal scalings 
are equivalent to the scaling of the standard quantum 
limit for unentangled atoms in the presence of transition 
broadening. 

The sensitivity of the continuous readout can be sur- 
passed if the system operates in a narrowband detection 
mode based on passive evolution, analogous to Ramsey 
spectroscopy. As first demonstrated in Ref. [El an equiv- 
alent passive Ramsey measurement technique can be re- 
alized using steady-state superradiance, interrupted for 
some period of time Td during which the decay rate 7 
and repumping rate W are set to zero. The atomic phase 
precesses during this dark period with no Schawlow- 
Townes phase diffusion and possibly lower systematic er- 
rors. Measurement of the light phase just before and after 
the shutoff allows the accumulated quantum phase to be 
estimated as (t>{Td) - 0(0) « VCTd) - V'(O) = The 
measurement record needed to estimate the final phase 
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tp{Td) also serves to estimate the phase for the next iter- 
ation of the experiment, potentiahy aUowing high repe- 
tition rate, non-destructive Ramsey-hke measurements. 

Ramsey spectroscopy detection operates with a detec- 
tion band centered at zero frequency. The favorable sen- 
sitivity of dc detection can be translated to in-principle 
arbitrary frequency using spin echo sequences to essen- 
tially serve the role of a mixer in a lock-in amplifiei^^. 
An even number of 7r-pulses must be applied in order 
to restore population inversion so that the superradiance 
resumes with the Bloch vector close to its initial steady- 
state inversion. 

Here we demonstrate the translation of enhanced sen- 
sitivity near dc to higher frequencies. Fig. [3^ shows a 
sequence in which the phase accumulation is coherently 
enhanced through a spin echo pulse sequence with the 
TT-pulses aligned in time to nodes of the applied modu- 
lation B{t). The phase difference Aip is approximately 
2/3n^, where /3 = 0.71 rad and rt^ = 2 for this data, 
and with only small corrections due to the finite 7r-pulse 
times of 15 fj,s. For a modulation phase that maximizes 
the accumulated phase deviation, the sensitivity (units of 
Tesla/VHz) in a 1 Hz bandwidth with no measurement 
dead time is 



V 



(4) 



where cr^ is the rms measurement noise of the light phase 
difference Aip. Fig. [Sj^ shows the phase advance Aip ver- 
sus the number of 7r-pulses applied for a 10 kHz modula- 
tion with /3 ~ 0.13 rad and modulation phase for which 
nodes of B{t) are aligned to the pulse times. The fit- 
ted slope is 0.28 rad/pulse, close to the expected slope of 
magnitude 2/3 = 0.26 rad/pulse. 

Fig. illustrates both that the sensitivity can be 
translated in frequency, and that the sensitivity increases 
with the number of 7r-pulses, while the bandwidth is 
decreased. This figure specifically shows the phase- 
insensitive transfer function 
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determined by measuring Aip and averaging over the 
modulation phase 9m at each modulation frequency fm- 
The spacing between the 7r-pulses is fixed to r = Td/nT^ = 
50 /iS, setting the maximum sensitivity to modulations 
near 10 kHz. Two different numbers of pulses are used: 
ri^r = 2 and 10 for the red and blue curves respectively. 
We find reasonable quantitative agreement between our 
data and modeling the effect of phase modulation by nu- 
merically integrating the Schrodinger equation for the 
two-level system in the presence of a modulated classical 
driving field with finite Rabi frequency D,/2t: — 33 kHz. 

In Fig.|4]3, we show that this technique allows sensitiv- 
ity below the limit imposed by PSN in the active mode. 
We compare the field sensitivity at frequency / = 36 kHz 
using an increasing number of spin echo pulses to the 
PSN-limited active sensing mode at the same frequency. 
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FIG. 4. (a) Normalized phase accumulation transfer functions 
measured for — {2, 10} and /3 = {0.25, 0.08} (blue squares 
and red circles, respectively), (b) Field sensitivity rj versus 
Td, calculated using Eq. [4] and the measured detection noise 
with average 0.45 rad. The sensitivity drops below the 
measured field sensitivity in the active mode (dashed line) at 
36 kHz. 



The phase measurement noise cr^ is roughly constant as 
Td increases, eventually leading to a sensitivity below the 
PSN-limited wideband value. The difference in scaling of 
sensitivity shows the gain in sensitivity in the Ramsey- 
like detection configuration at the cost of detection band- 
width. 

In conclusion, we have demonstrated a magnetometer 
that utilizes the active mapping of an atomic phase onto 
a light field phase to operate in both an active and passive 
field-sensing mode, with the choice of mode dynamically 
selectable. At low frequencies, the current magnetometer 
is limited to a minimum field sensitivity scaled to contin- 
uous operation of 190 pT/VHz. This is approximately 
30 times worse than the fundamental limit set by quan- 
tum phase diffusion, largely due to dispersive shifts of the 
cavity mode frequency, which in the future can be sup- 
pressed by operating at larger detunings not accessible in 
the current experiment, or by operation on two indepen- 
dent lasing modes for which the shifts are common-mode. 
Future experimental work will focus on moving to con- 
tinuous operation and realizing sensitivity at the phase 
diffusion limit. 
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